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ABSTRACT 


A pilot engaged in a strike mission against targets of opportunity 
is continually confronted with the decision as to how many and which of 
his available weapons should be expended against a variety of targets 
which are encountered. Associated with this decision problem is the 
problem of determining an optimal load for such a mission subject to 
the constraints of available payload and number of weapons stations on 
the aircraft. Certain assumptions are made concerning the distribution 
of targets within the target area which lead to a dynamic programming 
formulation of the decision problem. This yields a system of ordinary 
differential equations which are solvable recursively. In addition to 
the dynamic programming model, a sub-optimum determination of a "best 
load" is discussed. Although not as complete or precise as the dynamic 
programming method, this formulation is more readily adaptable to squa- 


dron level decision making. 
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1. INTRODUCTION 

The problem of selecting weapons for an aircraft strike mission 
and using those weapons effectively has been documented quite exten- 
sively. However, present publications are devoted almost exclusively 
to the situation where tle target is assumed known prior to the 
mission and the aircraft can be loaded for that target, or where 
several different aircraft can be loaded differently to cover an 


1 These approaches are inadequate for 


unknown spectrum of targets. 
determining weapons requirements against targets of opportunity. A 
pilot engaged in such a mission is confronted with the decision as 
to which and how many of his available weapons to expend when any 
target is encountered. When various types of targets might be en- 
countered, each with a different strategic or tactical value, he 
wishes to maximize the value of the targets he can expect to destroy. 
Associated with this is the problem involving the choice of weapons 
that should be initially loaded onto the aircraft for such a mission. 
Here the problem is constraíned by maximum payload, weapons, and wea- 
pons stations availability. It is desired to obtain, within these 
constraints, (1) the optimum feasible weapons mix to be loaded for 
use against an uncertain spectrum of targets, and (2) a doctrine for 
the expenditure of these weapons when targets are encountered. 

The model developed herein attacks these problems in reverse 
order. First, an optimum firing doctrine is generated for all feas- 


ible loads. Then that load which yields the maximum expected return, 


INAVWEPS Report 8664, "Preferred Aircraft Loadings Based on the 
Predicted Combat Effectiveness of Conventional Air-to-Surface Weapons,' 
NOTS TP 3710, 1965, (SECRET). 


when expended in accordance with the optimum doctrine, becomes the 


optimal load. 


ASSUMPTIONS IN THE MODELS 

Certain assumptions, necessary to the model development are made 
concerning target distributions within the target area and tactics of 
the strike aircraft. 

It is assumed that there are a finite number of different target 
types, each randomly and independently distributed within the target 
area with known density. For example, from photo-reconnaissance or 
other sources it might be known that the density of trucks within the 
target area is x per square mile. The number of targets of each type 
is assumed to be large enough so that the probability of future en- 
counters is unchanged by consideration of the number already met. 
Further, these targets are assumed to be of a unitary nature, charac- 
terized by the existence of single shot kill probabilities, as dis- 
tinguished from area targets characterized by probabilities of des- 
troying a certain percentage of the area or value of the target. 

The assumption is made that a ratio scale of target values can 
be specified. For example, with the least valuable target as a base 
with value V, the achievement of a kill on a target type i is worth 
7V, target type j is worth 8V, etc. 

For defensive reasons the strike aircraft is restricted to a 
single attack on any target which might be encountered. The pilot 
may not drop a bomb on the target, assess damage, and then make 
another attack if a kill has not been achieved. He may however, 


expend any number of weapons on this single pass provided they are all 
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the same type. That 15, while he may expend three 500 lb. bombs, he may 
not drop one 500 lb. bomb and one 1,000 lb. bomb. This restriction 
15 not necessary to the development of the models, but rather recog- 
nizes that differences in aim points between different weapons make 


this mixed drop undesirable. 


METHOD OF SOLUTION 

Using data on weapons weights, constraints of aircraft payload 
capacity, and the number of available weapons stations on the plane, 
all feasible loadings can be enumerated. Initial loads are then con- 
sidered which are not subsets Of some other feasible load. That is, 
if five weapons of type one and three weapons of type two constitute 
a feasible load, then four weapons of type One and three of type two 
is also feasible but a subset of the former. The kill potential of 
some subset of a feasible load is clearly less than that of the 
original load and is therefore not considered in the determination 
of an optimal initial loading. 

A functional equation for the expected value of targets killed 
is then set up as a dynamic programming process yielding a value for 
kill potential for each of the feasible loads being considered. The 
maximum kill potential from amongst these loadings determines the 
optimal load for the strike aircraft. Additionally in the solution 
of this dynamic programming process, the optimal firing strategy for 
the pilot is obtained for any initial load. The optimized sequential 
decision process yields a policy which determines as a function of 
the operating time and weapons remaining, the best weapons to be ex- 


pended in an encounter with any of the various target types. 
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Dynamic programming provides a complete and very precise means 
of determining the optimum load for a "targets of opportunity” type 
mission. Except in the very simplest cases however, the problem is 
quite complex and will require a computer for solution. In the majority 
of operational situations such equipment will not be availabie. Thus 
an approximate method, suitable for hand calculation, is needed for 
determining the initial load. 

In order to achieve simplicity, this "hand model" must sacrifice 
detail and thus completeness. In this respect it is a sub-optimal 
method. The simplified method is best suited for comparing a few, 
rather than all, feasible loads for use against a selected target 
list. By choosing the best weapon for each target type and comparing 
the expected returns for various mixes of these weapons, a "best load” 


can be generated. 
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2. GENERAL DISCUSSION AND DEVELOPMENT OF THE MODELS 


The random distribution of targets leads to an exponential inter- 
encounter time distribution for the strike aircraft searching the 
target reas (All targets which are within the visibility range of 
the aircraft's track are assumed sighted. Alternatively the sightings 
could be reduced by some ..,.eudation factor which takes into account 
the difficulty associated with perceiving a given target type. This 
however would not alter the analysis which follows.) The mean rate 
of sightings (the reciprocal of the mean of the exponential distri- 


bution of inter-encounter times) is denoted by r. 


DYNAMIC PROGRAMMING MODEL 
The development of the dynamic programming formulation of the 
problem is best seen by first considering the simpler cases. To begin 
with, a model has previously been developed for tle allocation of wea- 
pons over a limited time Т Е This can be utilized for gener- 
ating the optimal firing strategy in the simplest case of one target 
type-one weapon type. The optimal loading problem is not pertinent; 
for a single weapon type the best load is clearly the maximum number 


that can be carried on tiie aircraft. First define: 


f (t) = expected number of targets killed in the time 
n interval (0,t) if there are n weapons at time 
t = 0 and these are utilized optimally. 


1 
Appendix I 


2 

Bram, J., "Allocation of Weapons to Targets with Exponential 
Arrival Times in a Limited Time Intervai,' OEG Interim Research 
Memorandum 32, 1963. 
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An initial loading of N weapons and a time T allowable for search 
is assumed, and Е (Е) will Бе calculated for O < n < N and 0 < Е < T. 
Define: 


p = probability that j weapons expended against a target 
j will achieve a kill. 


Clearly £ (t) is dependent on the pes and on r, the encounter 
rate. 

By Bellman's "principle of optimality," for 0 < ts t, the 
optimal policy in the interval (0, t), given n weapons at time t = 0, 
must also be optimal in the subinterval (51, t) given n - n weapons 
at time t = ti» where 2^ weapons were expended in the first subin- 
terval (0, tj). 

Let E = At << T. By definition of the exponential arrival pro- 
cess, the probability of an encounter in a small interval of time At 
is proportional to At with constant of proportionality r, the en- 
counter rate. Thus the probability of one encounter in the time 
interval (0, At) is rAt and no encounters occur with probability 
ES TRAE 


Breaking the time interval (0, t) into the two sub-intervals 


(0, At) and (At, t) yields: 


(1 - råt) = the probability that no encounters will occur 
in the first sub-interval. 
p. = the expected return in the first sub-interval 


3 having expended j weapons. 


Ї 
Bellman, R., Dynamic Programming, (Princeton, New Jersey: 
Princeton University Press, 1957), p. 83. i 
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f (t - Ot) = the expected return in the second sub-interval 
having expended no weapons in the first and following 
an optimal strategy in (At, t). 


f | (t- Ot) = expected return in the second time interval having 
Пе] expended j of the n weapons in the first sub- 
interval and behaving optimally in (At, t). 


The expected return over the whole interval (0, t) is then: 
(1 - rôt) f (t - At) + rät (p, +£ (t - At)) 
n J n” J 


To optimize over the period (0, t), j must be chosen such that this 


expression is a maximum. Therefore 


(1) £ (t) = (1 - råt) £ (t = Aty + Tht max» (p FE . (t - ^t)) 
SSA i 


The maximization here is over the range lSj<sn rather than Osjsn 
since it clearly would never pay to let a target yo by without firing 
any weapons, saving them all for future encounters which are not 
certain to occur. 

Rearranging (1), dividing by At and taking the limit as At 


approaches zero yields: 


(2) f'(t) = - rf (t) +r max (p. +f (9) 
n n n-j 
lsjsn 

From (2), £ Ct) can be calculated for all n and t recursively 

as follows: By definition, f (t) = O for all t. (The expected 
0 

return is zero over any period of time if there are no weapons with 
which to operate.) Also f (9)= 0 for all n. (The expected return 


is zero if there is no time in which to operate.) Solving (2) for 
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п = 1 yields: 


БО) “р, (1-6 
Note that this result was to be expected.  Exponential inter-encounter 
times lead to a Poisson counting process for the number of encounters 
within a time interval. The term (1 - Ba), is simply the prob- 
ability of at least one encounter in the time interval (0, t), and Р, 
is the probability that, given an encounter, a kill will be achieved 


with the expenditure of the single weapon. 


Utilizing f (t) equation (2) can now be solved for n = 2. 


, зэ E rt 
2 = rf (t) + г max (p (2 e ), p) 


1 P | 
Since p > p. (2 - e!) when either t<- In огіғр > 2p,, 
2 1 2ру-р 2 1 
2 
the solution to the differential equation is: 
2 1 р 
Е (Е) =p, (l-e ES) τος t S ln IEEE. or if p 2 2p 
2 2 B 2p, -P 2 1 
12252 
-rt Р 
= Е Е + 4 - IT 
ЕЭ 2p, - e (rp,t + p, ΡΕ. σας ) 
172 
1 P1 
for t2-— Ш —— and p < 2p 
2р -р й : 
1222 


The significance in this solution lies in the "cross-over" point in 
time. If the search time remaining is less than this value of t, the 
optimal return is achieved by dropping both bombs on the next target 
which is encountered. If greater, only one bomb is dropped on the 


target, the other being saved for possible future encounters. 
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Having found f (t), 9 can now be calculated. Tlıe differential 
equation however becomes unduly complicated for n> 2. A numerical 
solution may be obtained for the problem though by recourse to equation 
(1), the difference equation, with which f (t) can be calculated for 
any n to any desired degree of accuracy by taking At small enough. 

With the possibility of encounters with targets of different 
types, the value of a target which is met must be considered when the 
decision is made as to how many weapons should be expended. 

The methods used ín developing the functional equation are 
analagous to those for the simpler model with the addition of con- 


siderations of target value and different kill probabilities. Define: 


m = the number of different target types. 


p = the probability that j weapons expended on a target 
type i will achieve a kill. 


r. * average encounter rate with type i targets. 
v = value of target type i. 


f (t) = expected value of targets killed in the interval (0, t) 
п given n bombs at time t = O and optimal utilization of 
these weapons. 


Then: 
m m 
(3) 01 See At) Е - At) ух At max (v,p,,tf  (t-^t)) 
f (t) " (1-, r, At) f, (t t LU | iP19*f. 
izl i=l Osjsn 


Note that in this case the maximization is over the range O < j < n 
vice 1 < j < п. Here the expected return might be greater if no 
weapons are dropped on targets of little value, expecting at some 


later time in the mission to meet something of greater wortn. 


The differential equation is obtained by subtraction of f (t) 
from both sides of the above equation, division by At and taking the 
limit as At approaches zero. 

m m 
£ (e) = | r, £,(t) +) x, ПАО сато? 


as] 1-1 0515п 


For n = l and m = 2 the solution to this equation can be shown tc be: 


(E P V t r,P5Vv)) Б 
J 1 1 222 (nee (r, + r.) t 


f (t) = ) for p v >p v 
1 (r. + r.) τον. 
1 2 
ап Ost< πιο ο. ας in 
-. Е p v = p. V 


or if pos 2 then for 


(r v ТҮ ν 
хээт M 1Р1Ў1 phos? 
E i 52 r (p v = p v ) 

12522121 2 2 


or if pv = Dos then for all t. 








l1 
L 
: Ё 
EDU ISD Vv DINE "r 
£ (t) =p v + (pv, еру) | 322 222 | сн 222 
1 11 2901 A E (p.v, - PV) ТЭР 
11 22 
r + 
1 111 ο. 
for d d and t Z Em 15 — r OU 
К? r. (P,V) ` P,V,)) 
21 
r F1P1Vi T r2P2V2 ] τη; “tt 
Е OA e 
2 92 P1" 
г.ртУүт Ж ГР. у 
for py Д and t 2 7 L ια T 
r - У 
1 2 ο. μι 
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Thus for the simple case of two target types and a single weapon, 
the expected returns can be calculated analytically. In addition, the 
optimal strategy is found, here consisting merely of a cross-over time, 
the time before which it does not pay to fire the weapon at the target 
with the lowest р Рог п > 1 however, it can readily be seen that 
solutions to the differential equation are not so easily obtained. 
Again though, through recourse to the difference equation, (3), numer- 
ical solutions can be obtaíned for any n to any desired degree of 
accuracy. 

Further generalization of the model entails a relaxation of the 
restriction to a single weapon type. The functional equation is 
Similar to those for the simpler models with some minor changes in 
the definitions of the terms. Define: 

р 2 the probability that k weapons of type j expended on 

ijk a target of type í vill achieve a kill. 
L = the number of different weapon types. 
f (t) = expected value of targets killed in the interval 


ΠΗ (0, t), given n, weapons of type i (for iml,---,L) 
at t = 0 and optimal utilization of these weapons. 


The difference equation can be written: 


m 
Е (O = (O E (t - At) 
non „on E d пп п ..п 
123 L 151 123 
m 
I АЕ max max бор... {ΕΤΕ} 
122 1 1 ijk an Th ee ell -k...n 
1-1 15151, “арал, 152 3 L 
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The differential equation can be written as before, but the 
solutions are obtainable only in the most trivial cases. 

A FORTRAN program has been written however from which numerical 
solutions to the difference equation can be σι οσα With inputs 
of the aircraft load carrying capabilities, the mission time, and the 


parameters of the model (m, E. v., L, p..,, and At) the program 
1 


ijk 

selects values for the n, which are consistent with the limitations 
i 

of the aircraft and maximize f. Additionally, another FORTRAN pro- 


gram! lists the optimal strategy to be employed in the decisions 


relating to the expenditure of these weapons. 


A SUB-OPTIMAL DETERMINATION OF LOAD 


In the majority of practical cases hand calculation will be the 
only means available for determining loading requirements. In such 
cases a sub-optimal method can be used at the squadron level for 
obtaining the initial weapons load for a mission. Define: 

T = search time (i.e., mission time allotted to searching for 

and attacking targets). 
р (Т) = probability that target type i will be met n times during 


ni the mission. 


со 


П (1) =}? ,9D = probability that target type i will be met at 
least once during the mission. 


n=1 
v = value of target type i. 
i 
m = number of target types. 
b, - number of sticks! of s weapons of type j in the load. 
5 


L = number of weapon types available. 


lin this development the terms stick and salvo are used equiva- 
lently. The choice of a stick or a salvo, in the load determination, 
will depend upon which gives the higher kill probability. 
20 


k = probability that a stick of s type w weapons will kill 
1W5 target type i. 


Note that in this development individual bombs are considered as 
a stick of one; i.e., two 250 lb. bombs might be considered as two 
sticks of one type j or one stick of two type j weapons. 

If, a priori, only the more highly valued targets and those 
individual weapons with the highest kill probabilities against these 
targets are considered] the assumption can be made that a target 
sighted will be dropped on, as long as the appropriate weapon for 
that target is available. (Note that this assumption is used only 
for determining the load. In the actual performance of the mission 
there is the possibility that, if that weapon most effective against 
a target is not available, a less effective weapon will be dropped.) 
The problem then becomes one of maximizing, for a given mission time, 
the expected value of targets killed, E, for those loads being con- 


sidered. E may be determined as follows: 
2 
E n f k V 
ντος ei 
1=1 


where £,. is the expected number of type i targets on which weapon 
ijs 


type j can be expended, and j indicates that weapon type for which 


у > 1 шоо eee . 
kiil E for all vw 123 ‚ 


If the number of type j weapons available were limitless, De 
со 


would simply be equal to a=) пр m. However, since a given load 
n=o 
is limited to b, sticks of type j weapon, no more than Pss can be 
js 


dropped regardless of how many type i targets are sighted. Thus 


тһе decision as to which targets and weapons will be considered 
will depend highly upon the judgement of the individual responsible for 
determining the aircraft load. 
21 


f... is given by 


ijs 
js = 
(4) O O 
.. =)n ) , 
ijs L Pai = σης 
n=o n=b 
js 


Equation (1) can be reduced to a more easily calculable form. 


Since _ 21 Бос! 
bis : 2 2 (D) ds 
| | Ὁ EON. ^N 
тмэис - 2 . ) Pat TE s 
п=о ™ n=1 n=b. n=1 
js 
equation (4) can be written as 
Т) om (+b D -b 2 Pay 
n=1 y ES JS ты] 
Б =1 
15 
τ 
=b H - (b. -n) p (Т) 
js i Bet 15 ni 


but II, 1 - p (T) and consequently b П -b -b,p. (T). Note 
= 0: js i js js Oi 


thatb p can be expressed as (b, - O) p . Thus equation (4) 
JS Oi Js Oi 
now reduces to b -1 
js 
B ale 
(5) =- (b, 7m) p (т) 
ijs js js ni 
n=o 


Equations (4) and (5) are valid for all b > 1. By definition 
35 


er = 0 for bjs = 0, since no sticks of type j weapon can be dropped 
if none are carried. 


As noted before, the random distribution of targets and conse- 


quent exponential inter-encounter times lead to a Poisson counting 
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process in determining tle probabilities of encountering given numbers 


of targets; i.e., for type i targets appearing at a mean rate r., 
i 


p .(T) = n EAD 
ni п! 


This gives b -1 
x Tm 
τ -г 
(6) f =b - (b,, 70) e i | 
n=o 


Table 11, Appendix III gives values of s for various com- 
binations of d and E This table greatly facilitates calculation 
of a "best load" from the feasible loads being considered. 

To determine the "best load' for a given set of targets, find 
that individual weapon type j, for each target, which has the highest 
kill probability. The loads considered are then the different mixes 
of these weapons selected from a list of feasible loads for the air- 
craft. 

For a given load, Boe will be known for each type weapon in the 
load and the encounter rate E will be estimated for each of the 
critical targets. The planned search time T will also be known for 


the mission. Thus using known values of τ. and estimated values 


of UM E can be calculated. 


Consider first the case where weapons will be dropped indivi- 
dually (sticks of one only) and where each weapon type in the load 
will be used against only one type target. In other words, target 
types and weapon types are matched one to one. Since weapons are to 
be dropped singly, Pe = P for each weapon type j. To determine E 
for such a load, first determine p and r.T for each weapon-target 


J 
type combination. Then obtain m from Table 11, and multiply this 
1] 
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value by k. v. Let £. k v =E . Add the resulting E 's to 


ijs i πο οσα ij Íj 
m 
get E for the load, 1.е., Е =\ g..- (Note that the E..'s are summed 
22503 1] 
i=l 


over i only, since the subscript j simply indicates that weapon which 
is most effective against target type i.) In the same manner, E may 

be determined for each load being considered; that load which yields 

the maximum E is then the "best load." 

For thís sub-optimal model, the extension of the problem to 
include sticks of two or more weapons is accomplished quite easily, 
since different size salvos or sticks of a given type weapon con- 
stitute different loads. For example, as noted previously, two 250 
lb. bombs might be considered as a load consisting of two sticks of 
one bomb or as another load consisting of one stick of two bombs. 
With this in mind, E can be calculated for these loads ir the same 
manner as previously outlined. In order to keep the problem hand 
workable, this approach assumes that each delivery of a given type 
weapon will consist of the same size stick. This assumption is made 
solely for determination of the aircraft loading prior to flying the 
mission. During the actual mission, however, it is quite likely 
that sticks of different sizes will be expended at different times. 
For example, a pilot might drop one type j bomb on the first truck 
he sees, a stick of three type j bombs on the second truck he sees, 
etc. Prior to the mission, however, this load of four type i weapons 
could have been evaluated as one stick of four, two sticks of two, or 
possibly as four sticks of one, since calculations for sticks of 


varying sizes quickly exceed the bounds of practical hand workability. 
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То addition to the case where targets and weapons are matched one 
to one, consideration must be given to the case where a given weapon 
type is optimal for more than one target type. ТЕ the assumption that 
a target sighted will be dropped on were adhered to, the expected 
nunber of type i targets killed would be proportional to the encounter 
rate m However, in this situation, a concession is made to the fact 
that, early in the mission, a target of lesser value might be passed 
in hopes of encountering a more highly valued target later. Thus it 
is assumed that the expected number of type i targets killed is pro- 
portional to both the encounter rate, Κι, and the target value, E 
The expected return from a given load is now calculated exactly as 
before with tl:e following minor modification: For all targets type i 
against which weapon type j is most effective, calculate the Ei 5 
exactly as in the simpler case (as if all type j weapons were to be 
expended on each type target). Then to approximate the particular 
expected value of type y targets killed by weapon type j, Бол 15 ши1- 
tiplied by the normalizing factor > = у where the summation 

3 riv 
in the denominator is taken over those type i targets for which weapon 
type j is optimal. These expected values are then added to give E as 
before. 

A logical extension to this hand approximation is the consider- 


' weapons. A secondary weapon would be considered 


ation of 'secondary!' 
to be a weapon which has high kill probabilities against several 
targets, yet is optimal against none of them. There are certainly 
cases where such a weapon might give higher expected returns over the 
entire target spectrum than would a load consisting of only those 


weapons which are optimal for the individual targets. For example, 
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consider a case in which there are two targets and three weapons to 
be used against these targets. The single stick kill probabilities 


for the individual weapons are as follows: 


WEAPON TARGET 1 TARGET 2 


1 «09 .01 
2 .01 .09 
3 . 08 .08 


If the load were limited to those weapons which are optimal for 
individual targets, only weapons l and 2 would be considered. However, 
it is obvious that weapon 3 merits examination, and, in fact, it appears 
as if a pure load of weapon type 3 might possibly be the best load to 
use in this particular case. Calculation of E for such a secondary 
weapon is identical to that in the foregoing case in which one weapon 
is considered for use against more than one target. 

This sub-optimal development may, at first glance, appear to be 
long and cumbersome. It must be realized though that in most actual 
situations all but a very few feasible loads will be eliminated from 
consideration prior to any calculation, and these will be examined 
only as far as they pertain to a small number of critical target tyres 
(normally not more than three or four). Thus it can be seen that 
relatively few calculations will actually be required and these will 
be of a trivial nature. 

In order to clarify the use of the sub-optimal model, the follow- 
ing simple example is presented in detail. 

Consider an armed reconnaissance mission with a planned search 


time of one hour. (Thus T = 1 so that rT ZI for all targets). 
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Three target types are considered: 


= 2 - 25 

у = ] r = 3 
2 2 

у = 4 г = 2 
3 3 


The aircraft has two weapon stations and two types of weapons 
are considered, one more effective against target type l, the other 


more effective against types 2 and 3 (see Table 2). 


The possible loads are indicated in Table 1. 










| WEAPON 


1 2 | 
Load | Single Stick of 2 | Single Stick of 2 





I 
| | 
| В 2 
| c 2 
| | 
D | 1 | 
Е 1 





TABLE 1 
COMPOSITION OF POSSIBLE LOADS 








* .06 1.08 | .03 
———————— 
1 01 1.05 1.02 
2) | 
2 .03 №12 1.051 
TABLE 2 


KILL PROBABILITIES 
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Load A: b =1,b =1 
11 21 


Target l - Меароп 1 


f = .99326 (from Table ll; b, - l, r T = 5) 
111 JS i 
E =£ К т = (.99326) (.02) (2) = .0397 
11 111 111 1 
Target 2 - Weapon 2 
£ = .95021 (b = 1, r T = 3) 
221 js i 
E =f К у = (.95021) (.05) (1) = .0475 
22 221 221 2 
Target 3 - Weapon 9 
£ = ,096&66 (b = 1, r T = 2 
321 JS 1 
= = (.86466) (.02) (4) = .0692 


E a: 
32 321 321 3 


Since weapon type 2 is optimal for two target types, E and E 


22 32 
must be normalized. 
NE «(у/су *rv)E = (3/11) (.0475) = .0129 
N - (8/11 NU = ,0608 
332 (8/11) (.0692) 


Thus for load A, E = .0397 + .0129 + .0608 = .1134 


Load B: D = 2 
11 
£ = 1.95283 (b. = 2, r. T = 5) 
111 JS i 
E = E (Load B) = (1.95283) (.02) (2) = .0724 
Load C: b = 2 
21 
£ = 1.75106 (b. = 2 r T = 3) 
221 JS i 


28 


NE, = (3/11) (1.75106) (.05) (1) = .0240 
f = 1.45866 (b. = 2, z = 2) 

321 JS 1 
НЕ = (8/11) (1.45866) (.02) (4) = .0846 
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E(Load C) = .0240 + .0846 = .1086 


Load D: b = 1 (Note that this is almost the same as Load B, but 
12 


now the two weapons are considered as one stick of two vice two sticks 


of one). 


f = .99326 
112 


Ed E(Load D) -» (.99326) (.06) (2) = .1192 





Load E: b = 1 


22 
f = .95021 

222 
МЕ, = (3/11) (.95021) (.12) (1) = .0312 
£ = .86466 

322 
NES, = (8/11) (.86466) (.05) (4) =.1252 


E (Load E) = .0312 + .1252 = .1564 

For the loads considered, the maximum E is .1564. Therefore the 
"best load" is load E which consists of one stick of two type 2 weapons. 
Note that if sticks of one only (case l in the model development) had 
been considered, the "best load" would have been load A with one of 
each type weapon yielding E = «1134. 

The foregoing sub-optimal method is limited in that it considers 


search time as a whole rather than a time-stepped firing doctrine as 


in the dynamic programming method. Thus those added values which 
might be generated by strategies which vary with mission time are 

lost in this formulation. Consequently the expected kill values 

will generally be much lower than those of the dynamic program formu- 
lation, and the "best load" generated will not always be the optimum 
load. However, the load obtained by this method should, if expended 
in accordance with the optimum firing doctrine of the dynamic program, 


be one of the more highly valued loads. 
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SAMPLE RESULTS 


3. 


OPTIMAL LOAD AND FIRING STRATEGY DETERMINATION 


For a mission time of 15 minutes, seven weapons from which to 


choose, and values of the other input parameters as listed in Table 


3, an optimal loading has been determined using the "BESTLD" program 


of Appendix II, with the following inputs: 


Weapon Type 


300 510 600 675 950 1975 
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Weapon Weight: 


Number of racks = 5 


Max payload = 5,000 lbs. 


Number of target types = 4 


Mission time = 0.25 hours 


Delta T = .OlO hours (This is the At of equation (3)). 
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KILL PROBABILITIES 
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These weapons' weights and aircraft payload and rack capacities 

result in approximately 400 possible loadings which must be considered. 
From amongst these 400, Table 4 lists the five loads with the greatest 
expected returns and for comparison the loading yielding the least. 

A loading designated 1200002 denotes one weapon of type one, two 
weapons ОЁ type two, zero weapons of types three through six, and two 


weapons of type seven. 


LOAD NO. LOADING EXPECTED RETURN 
1 0300002 1.29514 
2 1200002 1.25048 
3 0100102 1.22353 
ч 1000102 1.20670 
5 3000002 1.20409 
381 2011010 ‚38069 
TABLE 4 


LOADINGS AND EXPECTED RETURNS, 
DYNAMIC PROGRAMMING METHOD 


The expected returns from the optimal load are over three times 
higher than those of the worst. In fact the twentieth best load, 
0020300, had expected returns of only .77700. Thus the optimal loading 
yields returns more than 657, better than a loading which is better than 
95% of the possible loads. 

The above inputs and constraints have also been used in the sub- 
optimal model to generate a "best load." Weapon 7 is the best single 
weapon against targets l and 2, and weapon 5 is best against targets 
3 and 4. Thus the only loads considered in this formulation were 
0000500, 0000401, and 0000102. Table 5 lists six possible load con- 
siderations (including stick options) and the expected returns from each. 


32 


Loading E 


0000102 (Weapon 7 used as a stick of 2) . 5766 

0000401 (Weapon 5 used as a stick of 4) . 4458 

0000102 (Weapon 7 used as 2 sticks of 1) . 3824 

0000500 (Weapon 5 used as a stick of 5) . 3740 

0000401 (Weapon 5 used as 4 sticks of 1) . 3320 

0000500 (Weapon 5 used as 5 sticks of 1) ‚1432 
TABLE 5 


POSSIBLE LOADS AND EXPECTED 
RETURNS, SUB-OPTIMAL METHOD 


The "best load" using the sub-optimal hand method is 0000102 
considering weapon 7 as One stick of two. Although only the principal 
weapon for each target has been considered, note that weight constraint 
still allows the loading of one more of the lighter, less effective 
weapons (weapon 2 rather than weapon 1 since weapon 2 is dominant over 
the entire target spectrum). The load thus generated corresponds to 
the third most valuable dynamic programming load, and, if used in 
accordance with the dynamic programming firing policy, it would be 
expected to yield a return of greater than 907, of that given by the 
actual optimum load. 

With the best loading as generated by either method, the optimal 
strategy can be listed using the "POLICY" program of Appendix II. The 
results of this determination for loading 0300002 are shown in Figure 
l along with a format for displaying the strategy compactly. The 
time axis is mission time remaining and entries such 2x7 denote that 
two weapons of type seven should be expended. 

As an example, as the mission begins, the subloading is 0300002. 


If the first target which is met is a type one, and is encountered with 
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TARGET 
1 


2 


pov up 


Pa 











SUBLOADING 0300002 


г | 
Е 13287 J x 2 NO FIRE 
| 
| g'x 7 

ΣΕ 
р 
х $ x 2 
—T- 
O .05 10 2.1» 
SUBLOADING 0300001 
1х7 NO FIRE 
μμ 
о 
> x 2 d xX 
-一 一 一 一 一 一 | 

O .05 5110 „15 
ы eae: SUBLOADING 0500000 
EE хо? NO FIRE 

> 
3x2 
3 x 2 2 x 2 1 2 
ο .05 . 10 4. 5 
SUBLOADING 0200002 

2 

ов ЕР 1x2 NO FIRE 
2K 7 
а 2 

2x7 | 2 x 2 

O «052 ο 10 „15 
FIGURE 1 


OPTIMAL STRATEGY DISPLAY 
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«20 025 
«50 «25 


TARGET SUBLOADING 0200001 


Jt 1х7 NO FIRE | 


e 1 x 7 | 


2 | 2 x 2 
f BEEN NO 
2x 2 FIRE 
O .05 . 10 . 15 ‚20 .25 


BLOADINE 0200000 ___ ~ 


2 








2 2 xX 2 

5 2х2 

4 2 x2 1х2 №О ЕТВЕ 
— m с ο л аг 
ο . 05 „10 ame „20 „25 


SUBLOADING | 0100002 


а= МЕРЕ | 
ims | 2x NO FIRE 





1 

2 QNSE 

5 2x 7 

+ 2x7 | 175032 
0 .05 , lo .15 . 20 .25 
Т SUBLOADING _ 010001 

1 x NO FIRE 

2 1х7 


IN 
pa 
s 
N 


1x7 puso 
O «05 «10 525 20 „25 
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TARGET 


2. 


2 


1x2 


SUBLOADING 0100000 


NO FIRE 


1 x 2 


O . 05 


2x7 
Лр 


ο »02 


1х 7 


J Se 7 


ты 


S Pu ο, 
О . 05 


(The abscissa in these 





e 10 ol ο20 , 25 
SUBLOADING 0000002 
NO FIRE 
2 х 7 
2 х 7 


NO FIRE 


.10 .15 


SUBLOADING 0000001 


„20 «25 


NO FIRE 
ПЯ 


1 х 7 


NO FIRE 


| c NN 
. 10 E „20 „25 


figures is the time remaining in which to 


search for and attack targets.) 
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.22 hours remaining, the optimal strategy indicates that no weapons 
should be fired. If the next target which is encountered is a type 
four at time .17, the subloading still being 0300002, a drop of three 
weapons Of type two is called for. The sub-loading now is 0000002. 
If a type three target is now met, both of the remaining type seven 
weapons should be dropped. 

It is of interest to note the changes which occur over time in 
the optimal weapon to drop when a target is met. With the loading 
0300002 and target type one, with more than .08 hours remaining, no 
weapons should be expended. With .045 to .08 hours remaining, the 
optimal strategy calls for tne expenditure of one weapon of type two. 
Weapon type seven however has a higher kill probabilíty against type 
one targets. The type seven weapons at this point are being saved 
for possible future encounters with targets of type two and four, both 
of which are more valuable and occur frequently enough so that there 
is a significant probability of an encounter in the time remaining. 
Although weapon type two is the most potent against type three targets 
-the most valuable target type- the probability of meeting a type 
three target is rather small because of the low encounter rate of 0.10. 
It is for this reason that type two weapons are not being conserved. 
Later in the mission however, the optimal weapons to drop on a type 
One target are type sevens. At this point in time the probability of 
future encounters with more valuable targets is remote. 

Since the ultimate objective is tne maximization of the expected 
returns for the míssion, a consideration might be made for lengthening 
the search time at the expense of lower payload and weapon station 


availability. That is, perhaps an additional tank of fuel can be 
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carried on one of the aircraft's weapon racks, thus reducing the 
number available on which ordnance can be carried and also reducing 
the weapons payload, but increasing the search time. 

Table 6 lists the expected returns for the five best loads 
generated when the search time has been increased to 0.50 hours, the 
payload decreased to 4000 lbs., the number of weapon racks decreased 


to four, and the other input parameters unchanged from those of Table 


2: 
LOADING EXPECTED RETURNS 
0000002 1.78803 
0060301 1.42184 
0100201 1.30025 
0012221 1.29748 
0000400 1.28808 
TABLE 6 


BEST LOADS AND EXPECTED RETURNS 


Since the expected value of targets killed for the longer mission 
is significantly better than that for the shorter, 1.78803 compared to 
1.29514, it is clearly preferrable to make the trade-off of gas for 
weapons. 

An interesting feature of the kill probability matrix of Table 
3 is the preferred weapons with which to attack the most valuable 
targets. Weapon type seven is the most potent against type two tar- 
gets while weapon type five is preferred against type three targets. 
Target type three, althovgi the most valuable, has a very low encounter 
rate. This results i: tie fact that weapon type five does not appear 
frequently in the loadings which yield high expected returns for short 
missions. Letting the maximum payload equal 4000 lbs. and the number 


of weapons stations equal four, the preferred loadings of these two 
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weapons alone can be examined as a function of search time. The three 
loadings which must be considered are 02, 31, and 40, where now, 02 
denotes a load consisting of zero weapons of type five and two weapons 
of type seven. Calculations show that an initial loading of 02 is 
optimal for mission times of less than 1.6 hours. For mission times . 
of 1.6 to 8.9 hours, loading 31 is optimal, and finally for T > 8.9 


hours, loading 40 is preferred. 


COMPARISON OF EXPECTED RETURNS FROM NON-OPTIMAL 
FIRING STRATEGIES 


The gains in kill potential realizable through the utilization 
of the optimal firing doctrine can be demonstrated by comparisons with 
the expected returns from other strategies. As an example, consider 
the following case witlı non-optimal strategies: 
Let: 
Number of weapons = l 
Operating time = l hour 


No. of target types = 2 


Strategy A: Expend the weapon on the first target which is 
encountered. 


E, = the expected returns from strategy A. 

Strategy B: Expend weapon only if a target of type 2 is 
encountered. (Note that target type 2 is the 
more valuable). 


E, = the expected returns from strategy B. 


T.: The random variable describing the time to first encounter 
with target type i. (T, is exponentially distributed with 
mean l/r.) Е 
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= : < mi t, Ih + Prob. (T < min (t_, l hour 
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For several values of the parameters Р, and r , Table 7 shows the 
i 
results of calculations of the expected returns,using either of these 


strategies or the optimal. 


Expected Value 
of Targets Killed 





4 Gain Optimal 





| Strat. | Strat. | Optimal | Over 
| | А _ В Strat. | А В 
[р 3 
1 | | 
| | 
E | | 
| .963 |1.011 | 1.270 | 32% 267, 
r = 10 | 
1 i | 
E = 1 | | | ! 
2 | | 
| = ,3 | | | 
Р = .4 | | 
ux ο! 6 | .35%4 [1.018 | 1m | 185% | 
r = | | 
paa С. | 
2 | | | | 
р = .1} | | | 
| 1 | | | 
р = „4. | | | 
| 2 | 414 | 1.011 | 1.064 1507 | 5 | 
rj 710 | | | | | | 
| | | | | | 
га?! | | | 
2 | | | | | 
TABLE 7 


GAINS IN EXPECTED VALUE OF TARGETS KILLED ACHIEVED BY 
UTILIZING THE OPTIMAL FIRING STRATEGY 
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These results are representative of those obtainable using other 
values for tle input parameters. Generally, when one of these non- 
optimal strategies approaches th» optimal in expected returns, the 
other yields very poor results. When the non-optimal strategies yield 
roughly equivalent returns, the optimal dominates both by a significant 
amount; this is demonstrated in the first example of the table by gains 
of 26% und 32% through use of the optimal strategy. 

Comparisons can be made in this way with any conceivuble strategy. 
It seems reasonable too that techniques which consider the target 
values and encounter rates could yield returns comparable to those 
achievable through use of the optimal strategy. The determination of 
these strategies however would probably involve calculations no less 


complex than those of the dynamic program. 


SENSITIVITY OF EXPECTED RETURNS TO TARGET VALUES 


A complete parametric analysis is beyond the scope of this paper; 
however, some further comment must be made about target values. The 
developments of both the dynamic progrumning model and the hand appro- 
ximation are based upon tlie assumption that a ratio scale of values can 
be assigned to all possible targets. Such a determination is completely 
subjective and must rely solely on the judgement of experienced planners. 
Quantitative target values might be determined by presenting a list of 
targets to several military commanders and asking them to rate these 
targets according to military value. Targets could be assigned numbers 
consistent with their importance relative to a target of least value. 
The results of such a large number of opinions could then be analyzed 


by methods ot sensory psychophysics in order to assign representative 
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values to all targets. Methods for accomplishing such an analysis 
have been reviewed and consolidated by S. S. Stevens. 

The above method will obviously generate strategic values. These 
will certainly not be valid for every tactical situation (e.g., a tank 
about to overrun friendly troops is more important than a railroad 
bridge at that moment even though the bridge may have a higher strategic 
value). A strategic scale will, however, provide a starting point -- 
a base which can be adjusted to meet specific tactical situations. 

An analysis of the effects on the expected returns of incorrect 
selections of target values can be made rather easily for the simple 
case Of a single weapon and two target types. Assume the "true" ratio 
of the target values is a Now suppose that an erroneous ratio, 


V., has been used to calculate the strategy. The expected returns 


г? 

which result from the utilization of this strategy - these returns 

being computed using the true value ratio, V * - can be compared with 
Ї 


the maximum of the expected returns, this maximum being obtained when 


the strategy derived using V * is utilized. Define: 
r 
V 7 v/v 
r 2 1 
and let the mission time = 1 hour. With a single weapon, a strategy 
consists merely of a cross-over time, the time after which the least 


valuable target should be attacked if met. These times can be calcu- 


lated from: 
С]ру + гәр), 








t = — ноос) i£v »1 
1 2 275208 1 

шин = ie ET ifv <1 
r] To r (Pi = p2V,) 

t =o if Ve = 1 


15.5. Stevens, "A Metric for Social Concensus," Science, vol. 151, 
4 February 1966, pp. 530-541. 
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Using these cross-over times then, the expected returns, E, are cal- 


culated from: 




















О F i - (rtr) Ρ1Σ1 | -(r,+r,) sl 
Р, r (1 - e ) + 21 (1 -е ) t Our 
1 2 1222 
-τ (1-20) -r4(1-t) -(r.+r)t r 
2 
в=ру* (1-е ? + e (1 - e 12) < ) 
2 工 шш 
-r (l-t -(τ.Ε t r 
tpe 7 Улс 1 г) ) : if t < 1 hour, and V > 1 
1 r +r r 
1 2 
oon ur D. ΕΝ `(rytr2)t τι 
din € °: ? rytr, 
-r.(1-t) -(r +r )t r 
+руже (1-е l Z у—2— if t <1 hour, and V_ <1 
2 г r. +r r 
1 2 
Letting: 
р "p = .10 
1 2 
r * 2/hour 
1 


r = l/hour 
2 


Table 8 lists the results for several possible values of V *. 
r 


While the scope of these investigations of value sensitivity is 
obviously limited, it appears that Over-estimation of the value of 
targets of higher worth might impose less risk than under-estimation. 
Further examination in the area of sensitivity where several weapon 
types are employed simultaneously would be necessary before conclusive 


results could be stated. 
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ТЕ У is misestimated by a factor of z, 


the expected returns decrease by y. 


js 2=10 2=5 2=1/2 2=1/10 
| y y y y 

| 1 287, 22%, 5% 8% 
| 2 x 87, 67, 107, 377, 
(3 | ow | o |o | sm 
| 5 | 27, | 1% 3% 67% 

| 7 17, | 1% 17, 66% 
E 1% | 22290001: | 407, 

| TABLE 8 


PERCENTAGE DECREASE IN EXPECTED RETURNS 
DUE TO MISESTIMATION OF Vee 


++ 


4. CONCLUSIONS AND RECOMMENDATIONS 


A method for determining an optimal initial loading for missions 
against targets of opportunity has been developed. As a by-product 
of this load determination, an optimal strategy for the sequential 
expenditure of weapons has been found. The gains in the expected 
value of targets killed, when the initial loading is optimal and when 
an optimal strategy is employed, have been shown to be significant 
when compared with initial loadings and firing doctrines derived by 
the application of simple intuition and judgement to kill probabilities, 
target values, and encounter rates. 

The determination of an optimal loading and expenditure policy, 
however, are at present of such complexity as to require the use of a 
high speed computer in all but the most trivial cases. For those 
operational situations in which a computer will not be available, a 
hand-calculable approximation method for the determination of a best 
initial loading has also been developed, which, though sub-optimal, 
has been found to generate loads which yield high expected returns. 

Either of the above methods should yield significantly better 
results than present policies which prescribe that aircraft be loaded 
for targets which are assumed known or that several aircraft each be 
given different loads in order to cover an unknown spectrum of targets. 

The models presented in this paper assume that weapons can be 
loaded in any combination which satisfies rack and weight constraints, 
and that the pilot has complete freedom of choice as to the sequence 
in which weapons can be expended. For most attack aircraft these 


assumptions are not always true. Restrictions on unsymmetrical loadings 
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are one limiting factor. Others relate to the mechanics of multiple 
bomb racks. The models should be expanded to account for these aspects 
of the problem. 

A vital assumption in the developments of both models was that 
inter-encounter times would be distributed exponentially. It is 
recommended that an analysis of operational data be made to ascertain 
the validity of this supposition. 

It seems reasonable to assume that a list of all possible types 
of targets of opportunity, for a spectrum of scenarios, would be 
small enough - certainly not larger than 50 - that a thorough analysis, 
as suggested in Section 3, could be made of their relative values. Such 
an analysis would indicate general strategic values which could then be 
adjusted for specific tactical scenarios. Furthermore, additional 
investigations similar to that outlined in Section 3 should be under- 
taken to study the sensitivity of the optimal load and strategy cal- 
culations to these values. 

The hand approximation method is rather crude in form. It is 
recommended that this model be examined further in an attempt to refine 
the calculations and to extend the model to cover a wider variety of 
cases in a more complete fashion. 

The computer programs of Appendix II were written solely for the 
computations used in this thesis and should be improved upon before 
extensive investigations are made of the dynamic programming treatment 
of the problem. The four target-seven weapon example required two 
hours and fifteen minutes on a CDC-1604 computer for the determination 
of an optimal load. This is testimony to the inefficiency of the 


"BESTLD'" FORTRAN program. 
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The problem of determining an optimal initial loading is similar 


l of dynamic programming. In the 


to the classic "flyaway-kit" problem 
strike aircraft loading problem the additional complication is present 
that the costs of shortages are not constant but depend on the initial 
loading. It would seem though that research into solution of the 
strike aircraft problem by this method should be fruitful. Very good 
approximations to an optimal initial loading should be obtainable very 
rapidly in this way. 

Although the "POLICY" FORTRAN program computes the optimal 
Strategy with relative ease, usually less than four minutes is required 
On a CDC-1604 computer, it is recommended that methods of approximation 


be developed for which a strategy can be calculated without the aid of 


a high speed computer. 


1 
Bellman, R. E. and Dreyfus, R. E., Applíed Dynamic Programming, 


(Princeton, New Jersey: Princeton University Press, 1962), pp. 42-43. 
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APPENDIX I 


ENCOUNTER RATES 

In a combat situation, many targets (especially "targets of 
opportunity") may be thought of as being distributed randomly through- 
out a given area, in that the pilot flying the mission does not know, 
beforehand, where or when a target will be encountered. This may be 
thought of as the classical situation in which it may be assumed that, 
in any small time interval dt, the probability that a target is en- 
countered is proportional to dt and is independent of any previous 
encounters. 

Over a period of time, photo or visual reconnaissance can give 
an average density of targets in a given area. From the attack air- 
craft speed and the average search sweep width (width of a path in 
which the pilot could be expected to see a target) the area swept per 
unit time can be obtained. By multiplying the average type i target 
density by the average sweep area per unit time, an estimate of the 
average encounter rate r. is obtained for target type i. 

The average encounter rate developed above can be used to gen- 
erate the distribution of inter-encounter times. 

Assumptions: 

Probability of one encounter in dt = r.dt + h(dt) where h(dt) 
approaches zero as dt approaches zero. 

Probability of two or more encounters in dt = h(dt) 

Probability of no encounters in dt = l-r,dt + h(dt) p(t)dt = 
prob. (t S first encounter time S t + dt) = dP(t} 


Е 
P(t) - | p(T)dt = prob.(encounter time < +) 
ο 
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Q(t)? = 1-P(t) 

Thus if t. is the encounter time, and t is any arbitrary time in 
the mission, 

Q(t+dt) = prob.((t > t) (no encounter in dt)) 

= Q(t) (1-r,dt + h(dt)) 

Subtracting Q(t) from both sides, dividing through by dt and taking 
the limit as dt approaches zero gives dQ(t) /dt = "Gr, 
Rearranging gives dQ(t)/Q(t) - d 
Integrate both sides to get ln Q(t) = Си 
Thus Q(t) = Me His 
But Q(0) » 1 since there can certainly be no encounters in zero search 


ет. 


: er 
time. Therefore A = l and Q(t) =e + and consequently P(t) =l-e i. 


and p(t) = reif 
Thus it ís seen that the ínter-encounter times are exponentially 
distributed. It also follows naturally that the probability distribution 
on the number of targets type i encountered during a mission of length 
T is в Lies, 
n -r T 
prob. (N(T) = n) = (r T) ETE the 


and the expected no. N(T) - r T (This is obtained from the observations). 
3 


за E., Modern Probability Theory and Its Applications, 
(New York: John Wiley and Sons, 1960) p. 253. 
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